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Abstract 



We have derived the general solution of a wave equation describing the 
dynamics of two-layer viscoelastic polymer materials of arbitrary thickness de- 
posited on solid (quartz) surfaces in a fluid environment. Within the Voight 
model of viscoelastic element, we calculate the acoustic response of the sys- 
tem to an applied shear stress, i.e. we find the shift of the quartz generator 
resonance frequency and of the dissipation factor, and show that it strongly 
depends on the viscous loading of the adsorbed layers and on the shear stor- 
age and loss moduli of the overlayers. These results can readily be applied 
to quartz crystal acoustical measurements of the viscoelasticity of polymers 
which conserve their shape under the shear deformations and do not flow , 
and layered structures such as protein films adsorbed from solution onto the 
surface of self-assembled monolayers. 
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1. INTRODUCTION 



The analysis of viscoelastic behavior of biomolecular films at fiuid-fiuid and fluid-solid 
interfaces is a very active interdisciplinary field and a "hot" topic in the physics of modern 
surface phenomena, which has many possible applications [1-13]. Particular attention is 
focused now on macromolecular multilayer assemblies such as polymer "brush" films [2], ad- 
sorbed protein layers [3,6,7-10] as well as on self-assembled monolayers (SAM) [1,6,9]. The 
importance of studying such biomolecular sandwich structures arises from the prospect of 
their use as molecular biosensors and for medical and enviromental applications [3,6,8-10]. 
The mechanical properties of these structures can be determined with high accuracy from the 
frequency response of a quartz oscillator in contact with them [6,8-14]. An increasing num- 
ber of theoretical and experimental investigations of quartz crystal microbalances (QCM:s) 
functioning in liquid enviroments have opened a possibility for in situ measurements of 
biomolecular assemblies and for a quantitive interpretation of their viscoelastic response [7- 
17]. Among theoretical methods, a continuum mechanics approach is often used. In contrast 
to the electrical circuit method or a transmission line analysis ( for a review see also 

[8] and Refs. [13,15] therein), this approach is directly linked to the physical description 
of the system and provides a clear connection between measured characteristics and ma- 
terial parameters. Within continuum mechanics the mechanical properties of viscoelastic 
materials are usually related to energy storage and dissipation processes resulting from the 
balance between an applied stress and the ensuing relaxation processes in the material. The 
experimental conditions determine the amount of stress, whereas the relaxation rates are 
intrinsic properties of the material [2,7,18,19]. 

The simplest way of accounting for the mechanical properties of a visco-elastic material 
is to introduce a shear viscosity coefficient rj and a shear elasticity modulus fi within one of 
two basic models due to (i) Maxwell and (ii) Voight. In contrast to the Maxwell model, the 
characteristics of the so called Voight element used in that model is the following: It does 
not describe a flow at a steady rate; the viscoelastic element is described by complex shear 
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modulus, real part of which (the storage modulus) is independent of frequency whereas its 
imaginary one (the loss modulus) increases linearly with frequency [19]. In addition to the 
Maxwell and Voight models various combinations of them have been used (for a review, see 
Refs [14,18,2,7]). 

The Maxwell model is usually applied to polymer solutions which can — at least for low 
shear rates — demonstrate purely liquid-like behavior. The Voight model is applicable for 
polymers which conserve their shape and do not flow [18,19]. 

In the present study we analysed, within the continuum mechanics approach, the vis- 
coelastic response of two polymer layers covering a quartz plate subject to shear defor- 
mations. The viscoelastic material was modelled as a Voight element, with its parallel 
arrangement of a spring and a dashpot. Even this relatively simple theoretical model has 
allowed us to demonstrate the most important consequencies of adsorbate viscoelasticity in 
terms of a deviation from the Sauerbrey relation [^ relevant for a gas phase environment 
and ultrathin layers. We have also been able to elucidate the role of viscous loading when 
the oscillator is immersed in a bulk liquid. 



We consider the case of two viscoelastic layers covering the surface of a piezo-electric 
plate oscillating in a pure shear mode in a bulk liquid. We derive here a general solution 
of the corresponding wave equation and analyze it using "no-slip" boundary conditions. In 
our model the result of Reed et al. [16] emerges as a special case (i.e., when both layers are 
equivalent, and the upper (Newtonian) liquid has a vanishing density). 

The geometry of the model system is shown in Fig. 1. We treat the quartz plate (index 
'0') as a harmonic oscillator whose resonant frequency in vacuum [ 4, 5, 10-17], /o, is 



Hence /o is a function of the thickness Hq of the quartz plate, its density po, and its elastic 
shear modulus /io- 



2. MODEL 
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The quartz slab is treated as a loss-less generator when it operates in vacuum (or a 
gaseous environment). It is known that mass loading changes the resonant frequency of the 
oscillator. For sufficiently thin non rigid overlayers, the shift of the resonant frequency is 
proportional to the added mass only [4] , 

A/ Ki foAm , 

and the energy dissipation due to viscous losses is then negligibly small [6] , i.e. such a layer 
demonstrates a solid-like response in vacuum or gaseous experimental conditions. 

In the opposite case of total immersion of the quartz slab into a liquid (in general into a 
viscoelastic medium) both the energy dissipation due to the oscillatory motion excited in the 
interfacial solid-liquid region, and the resonant frequency shift, are functions of the overlayer 
viscosity (elasticity) and density [13]. Our results show that even very thin viscoelastic 
films can dissipate a significant amount of energy when the quartz plate with overlayers 
oscillates in the liquid phase. The energy dissipation in QCM experiments can be obtained 
by measuring the dissipation factor D which is inversely proportional to the decay time 
constant [6] 

TT/r 

where / is the resonant frequency. For a purely viscous bulk liquid, it has been shown that 
the resonance frequency shift [15] and the dissipation factor shift [5,6] depend on the density 
and shear viscosity of the hquid. 

Let us consider a viscoelastic overlayer medium within the Voight scheme of using a 
viscoelastic element (Fig. 2) which consists of a spring and a dashpot in parallel. When a 
shear stress aa:y is apphed to a Voight element, the elastic response of the spring and the 
dashpot viscous resistance contribute to the stress/strain relation as follows: 

= + • 

Here Ux is the displacement in the x direction and the corresponding velocity is Vx] n is the 
elastic shear modulus and r] is the shear viscosity of the overlayer, respectively. 
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The wave equation for bulk shear waves propagating in a viscoelastic medium is [16]: 



(2) 



where = + ifj," — 11 + iurj is a complex shear modulus. The general solution of Eq. (2) 
can be written as 



u,{y,t) = {C,e-^y + C2e^y) e'^\ 



(3) 



where ^ — a + ik contains a decay constant a, 



a — 



(4) 



and a wave number k, 



s\ 1 + 



X = : 

■quj 




(5) 



Here we introduce the viscoelastic ratio x the ratio between real part (storage modulus) 
and imaginary one (loss modulus) of complex shear modulus and a viscous penetration depth 
5, respectively. The coefficients Ci and C2 in Eq. (3) can be determined from the appropriate 
boundary conditions. 

3. "NO-SLIP" CONDITIONS 



The so-called "no-slip" boundary condition at the solid-overlayer interface {y — 0) cor- 
responds to a continuous variation of the displacement Ux and of the shear stress across 
the interface. These quantities are related as follows, 



_ dqx _ 

oy 



(6) 
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Here Qx correspondes to the component of quartz substrate displacement vector, //q is the 
substrate elastic shear modulus. 

For the two-layer system we have a continuous variation at the surfaces y = hj of the 
velocities v^^ and the shear stresses cr-j.. Hence at y — hi we have 

(1) (2) 

v'^'=vf\ .i^) = .f = 0, 

4' - -P.5.k + + ^) + /^i I ^ + ^ 1 , (7) 

The no-slip boundary condition at y — h2 gives: 




(2) (3) 

Finally, at the free surface, y — hs, the boundary condition is 

5^,(3) 

A general solution for the wave equation (2) with boundary conditions (6-9) can be written 

as 

e2ay + ^e2^i/vi 

_ (l + ^e^^^^^^) - (l - ^e^^^^^-^) K:,6 + ^3e3tanh(e3A/^2) . . 

«i6 (1 + /e^^^'^O + «26 (1 - /e^^^'^i) ' - ^6 tanh {^s^h2) ^ ^ 

A/li = /l2-/il, A/l2 = /i3-^2, (12) 
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«.,2 



"1,2 + iki,2, 6 



( 

V3 



(13) 



«1,2 = ?7l,2 — = -J-^, «3 = ?73 



UJ 



(14) 



The shift in resonance frequency A/ and dissipation factor AD are obtained, respectively, 
from the imaginary and real parts of the /3- function [6] : 



A/ = Im ^ 
AD = -Re 



/5 



yT^fpohoJ 



(15) 



where 



/3 = Ki^i 



1 - Ae^^i'^i 



(16) 



1 + Ae^^i/^i ■ 

Below we analyze the general results expressed in Eqs. (10-16) in the hmit cases of thin 
and thick viscoelastic layers deposited on the quartz oscillator surface when the oscillator 
operates (i) in vacuum (or a gaseous environment) and (ii) in a bulk liquid. 



4. RESULTS 
4.1. Viscoelastic layer in vacuum 

Prom general solution for the wave equation (10) it follows that the overlayer thickness 
always appears as a dimensionless combination. Thus, the limit cases of " thin" /" thick" 
film correspond to the situations wheather the film thickness is much smaller/greater in 
comparison with the reverse value of a decay a constant and the propagation k constant 
which is reduced to the following criteria: 

in the 'thin film' approximation h/5 <^ ^/x/^ for x > 1 and h/S — x/2 for x < 1; 

'thick layer' correspondes to inequalities h/5 ^ X\/2x for x > 1 and S> 1 + x/2 for 
X < 1- 
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For pure viscous film this criterium corresponds to the well-known ratio between a film 
thickness h and a viscous penetration depth S: in the thin layer approximation h/S <^ 1 
[20]. 

Let us consider the solution of wave equation for one viscoelastic layer on the QCM 
surface when the system oscillates in vacuum (or in a gas). In such case, from general 
expressions (11-16) we get for /3-function: 

/?= -<tanh({/i). (17) 

By a series expansion on ha <^ 1 and hk ^ 1 in special case of thin overlayer, we find from 
(15) and (17) that 

^ - oJ!"! 2, ) (18) 

2TTpoho 3d^(l + x) 

AD ^ ^^'^^ ^ , ^=JL,S = ,i^ (19) 

It is clear from Eqs. (18) and (19) that to linear order in the (small) thickness, the 
dissipation factor vanishes, while the frequency shift is a function of oscillation frequency 
and mass of the layer. This is the well-known Sauerbrey relation for the loss-less quartz 
oscillator overlayered with a thin film: 

/o 

Sauerbrey = TT^'^^ = 

AD = 

The viscosity and elasticity of the overlayer only appears in the third order approximation 
with respect to thickness. For sufficiently thin layers, it leads to a small correction to the 
Sauerbrey relation which can be hardly observed experimentally. 

In the opposite case of a "thick" viscoelastic layer {ha 1, /i/c 1) we readily find 
that P —K,^ and that 
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(20) 
(21) 



The shift of resonance frequency (20) and the dissipation factor changes (21) are functions 
of film viscoelasticity and can be measured in QCM experiments. 

In the viscous hmit (// — > 0) we recover the result of Ref. [15] for the frequency shift, 

1 



A/ 



prjLj 



27rpo/io 

and the result of Refs. [5,6] for the shift of the dissipation factor, 

1 ' 



(22) 



AD 



'2pr/ 



(23) 



4.2. Viscoelastic layers in a bulk liquid 



Let us consider two thin viscoelastic ovcrlaycrs of thickness hj {j = 1, 2) under a bulk 
Newtonian liquid for which its thickness A/i2 satisfies the inequality ^3 A/i2 ^ 1 . In this case 
we may expand the expression for /3 and keep only terms linear in hjaj <S 1 and hjkj <S 1. 
The acoustic response of the QCM with the thin viscoelastic layers immersed in a Newtonian 
bulk liquid (index "3"): 

1 



A/ 



2-11 poho I ^3 
1 



E 

V3 



hjPjCU — 2hj 



V3 



AD ^ — = — { — + y 



2/V 1 1^ 



d's ) p,j+ u^r]^ 
2 



PjLU 



3/ p^j+LJ^f]] 



(24) 
(25) 



Prom Eqs. (24) and (25) it follows that for ultrathin films the contribution of the film is 
small in comparison with the bulk liquid acoustic response A/ ~ —ffs/S^. AD ~ tjs/Ss- 
However, a thin layer with a finite thickness will demonstrate a different acoustic response 
depending on the ratio between the viscosity and the elasticity of the film. 
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To make a picture more clear, let start from only one thin viscous overlayer under a 
bulk Newtonian liquid. For a purely viscous overlayer, the resonant frequency shift depends 
on a bulk T^a/^a term and the difference between the film mass contribution {hipiuj) and 
its viscous contribution (~ hi/rji) multiplied by factor {rj^/S^)^ due to the penetration of 
acoustical vibrations through the adsorbed layer into a bulk liquid and corresponding energy 
loss due to the viscous dissipation. The dissipation factor (25) is determined by properties 
of bulk liquid. 

We found, that even in this simple case there exists an important for experiments sit- 
uation when the resonance frequency shift of layered system (24) is equal to a bulk liquid 
responce only but not of adsorbed layer one. This occurs for such a relationship between 
viscosities and densities of the media, i.e. rjsps — rjipi, at which the 'mass responce' will 
compensate the viscous one. We can explain this subtraction as a result of damping and re- 
flection of travelling acoustic waves in a layered medium with different density and viscosity 
of the layers. 

In the opposite case, for purely elastic film (jji — > 0) one can find that only the dissipation 
factor change depends on the film elasticity (AD ~ /ii/ pi) , while the frequency shift contains 
the mass response of the film (A/ ~ /iiPi) but not the elastic one. It is easy to see that 
for a gaseous environment when 773 tends to zero, we recover the Sauerbrey relation (see 
also the next paragraph). Thus, it is essential to note that the viscous bulk loading of the 
adsorbed layer leads to the apparence of a viscous (viscoelastic) contribution of the layer 
and hence opens up a possibihty for the corresponding QCM measurements to determine 
the viscoelastic parameters of the overlayer(s) when the resonator oscillates in a bulk liquid. 

In Figs. 3 - 6 we have presented the calculated change in resonance frequency and in 
the dissipation factor as a function of the viscoelastic film thickness hi for a given value of 
the layer viscoelasticity {pi — 10'^ - 10^ N/m^, rji — 0.01 Ns/m?, pi — 10^ kg/m^) when 
a quartz plate oscillates with frequency /o = 5 10^ Hz in bulk water (p3 = 10^ kg/m^, 
7^3 = 10"^ kg/m^). 
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In real experimental situation, not only changes of the thickness of adsorbed layer must 
be taken into account. Variations of A/(/x, 77) and AD(/i, 77) during the adsorption can be 
of crucial importance as viscoelastic parameters of the adsorbate can also vary during the 
process of coverage. 

To illustrate the variation of A/ and /S.D with the viscoelastic layer parameters ji and 
77, we have shown in Figs. 7-8 the results of a numerical simulation of the model for the 
thin viscoelastic layer of a given thickness hi — 10~^m and density pi — 10^ kg/m^, when a 
generator oscillates with frequency /o — 10^ Hz in bulk water. The resulting nonmonotonic 
behavior of Af{ij,, rj) and AD{iJ,, rj) gives the evidence for key contribution of viscoelasticity 
in complex polymer adsorbed layers dynamics in a liquid medium. 

5. CONCLUSIONS 

We have found a general solution of the wave equation describing the propagation and 
decay of shear waves in layered viscoelastic medium in contact with a quartz surface. This 
sandwich structure can be used to model the response of viscoelastic polymer and biomolec- 
ular films such as polymer "brush" films or proteins adsorbed from solution onto a solid 
substrate modified by SAM covering before adsorption. Our results can be summarized as 
follows. 

(i) We have analysed both "thin" and "thick" layer acoustical responces. These two limit 
cases correspond to wheather the film thickness is much smaller /greater in comparison with 
the viscous penetration depth 5 — ^j2j]Jpuj and the viscoelastic ratio x = IJ'/v'^- 

(ii) To elucidate the validity of the Sauerbrey "ideal mass thin layer" response, we have 
calculated the "thin" viscoelastic layer resonance frequency shift. We found that the "ideal 
mass layer" response corresponds to a purely elastic behavior of the Voight element, i.e. to an 
element with zero shear viscosity. Due to the viscosity of the adsorbed layer, the acoustical 
responce of the system becomes to diverge with an "ideal mass thin layer" behavior when 
the ratio of the film thickness h and the penetration depth 5 tends to 1/2 (for small values 
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of viscoelastic ratio x ^ 1). 

(iii) In case of a "thick" viscoelastic medium we found the shifts in resonance frequency 
and dissipation factor as functions of the ratio of storage and loss moduli. Using Eqs. (20) 
and (21) it is possible to calculate these moduli using data from the corresponding acoustic 
experiments. 

(iv) We derived the general expression for the acoustic response of two viscoelastic layers 
of arbitrary thickness covering the surface of a quartz oscillator when it is immersed in a 
Newtonian liquid. This result allows us to analyze the dynamics of a layered medium with 
arbitrary layer thicknesses and with a viscoelastic ratio which can vary during the coverage. 

We have presented analytical expressions for the resonance frequency shift and the dis- 
sipation factor response of thin viscoelastic layers in contact with a liquid. The analysis 
showed that in contrast to oscillations in vacuum or in a gas environment, even a thin vis- 
coelastic layer will dissipate a significant amount of energy. We have therefore also shown 
that a viscous bulk loading of the overlayer opens up the possibility to measure its viscoelas- 
tic properties from simultaneous measurements of A/ and AD using a QCM operating in a 
liquid enviroment. 

The results of the model can readily be applied to quartz crystal acoustic measurements 
of adsorbed proteins which conserve their shape during adsorption and do not flow under 
shear deformation and of the polymer films far from the glass transition region. 
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FIG. 1. Geometry of a quartz crystal microbalance (QCM) covered by a double-layer vis- 
coelastic film. The QCM system oscillates in a bulk liquid. 
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FIG. 2. Schematic depiction of the Voight viscoelastic element. 
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FIG. 3. Numerically calculated resonance frequency shift as a function of a viscoelastic over- 
layer thickness hi; pi = 1000 kg/m^, rji = 10"^ Ns/m^, = 10^ N/m^. The QCM osciUates in 
bulk water; pi = 10^ kg/m^, r?3 = lO^^ Ns/m^, /o = 5 • 10^ Hz. 
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FIG. 4. Numerically calculated dissipation factor change as a function of a viscoelastic 
layer thickness (for the same parameters as in Fig. 3). The QCM oscillates in bulk water. 
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FIG. 5. Numerically calculated resonance frequency shift as a function of a viscoelastic over- 
layer thickness hi; all the overlayer parameters are the same as in Fig.3 but with /xi = 10^ N/m^. 
The QCM oscillates in bulk water; /o = 5 • 10^ Hz. 
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FIG. 6. Numerically calculated dissipation factor change as a function of a viscoelastic 
layer thickness (for the same parameters as in Fig. 5). The QCM oscillates in bulk water. 



20 



Shear modulus (Ns/m2) 



Viscosity (N/m2) 



FIG. 7. Numerically calculated resonance frequency shifts as a function of the shear elasticity 

modulus fii and the shear viscosity rji of the viscoelastic overlayer of thickness hi = 10~^ m. The 
QCM oscillates in bulk water; /o = 10^ Hz. 
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FIG. 8. Numerically calculated dissipation factor change as a function of the shear elasticity 

modulus fii and the shear viscosity rji of the viscoelastic overlayer of thickness hi = 10~^ m. The 
QCM oscillates in bulk water; /o = 10^ Hz. 
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